Local entanglement and quantum phase transition in ID transverse field Ising model 
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In this paper, we study the entanglement between two-neighboring sites and the rest of the 
system in a simple quantum phase transition of ID transverse field Ising model. We find that the 
entanglement shows interesting scaling and singular behavior around the critical point, and then 
can be use as a convenient marker for the transition point. 
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Entanglement is a one of the most fascinating features 
[H of the quantum theory. It roots in the superposition 
principle of quantum mechanics and the tensor product 
structure of different Hilbert spaces, and implies the exis- 
tence nonlocal correlation in quantum world. Since this 
correlation is absent in the classical world, it is a kind 
of pure quantum correlation and regarded as crucial re- 
source in many quantum information processing tasks 
0, 0] • Actually, similar fascinating property also exists in 
the classic optics. By recording the interference pattern 
of two laser beams reflected from one object, the whole 
image can be recovered from even a small piece of holo- 
graph, although the resolution may be reduced. In other 
words, the superposition of two coherent electromagnetic 
waves enables us to learn some global information from 
a localized spatial area. Though it is not a precise anal- 
ogy, people expect that the entanglement, which roots in 
the same superposition principle, can enable us to learn 
some global properties from a small part of the system. 
This observation may be one of the main motivations in 
the recent studies 0, i, 0, 0, & i, ES El El El on the 
role of entanglement between a small part, e.g. a block 
consisting of one or more sites, and the rest of the system 
in the quantum phase transition [Til ]. 

Conceptually, the single-site and the two-site entan- 
glement, and the block-block entanglement [ill, E3, EH 
are different notations of the same kind of measure- 
ment. They arc defined as the von Neumann entropy 
of the reduced density matrix of a local part in the sys- 
tem. Studies of the block-block entanglement in the one- 
dimensional spin models have established a close rela- 
tionship between the conformal field theory and quan- 
tum information theory [12]. The block-block entangle- 
ment is also generalized in systems other than spin mod- 
els. Many works are devoted to understand the local 
entanglement, a limiting case of the block-block entan- 
glement, in the ground state of itinerant fermion systems. 
For example, in the extended Hubbard model, the global 
phase diagram can be sketched out by the contour map of 
the single-site local entanglement in the parameter space 



0. And the entanglement's scaling behavior is investi- 
gated in the one-dimensional Hubbard model at critical- 
ity point [f|. In the Hirsch model, the singularity in the 
first derivative of the single-site entanglement is used to 
locate the superconductor- insulator transition point 0]. 
However, in some spin systems with Z2 symmetry, the 
single-site entanglement usually does not contain enough 
information to describe the critical behavior [1, Q , then it 
is necessary to include more sites, such as two sites in the 
simplest case, into the block. For example, in ID XXZ 
model, the two-site local entanglement reaches maximum 
at the isotropic transition point [9(. 

These results suggest that the local entanglement may 
be used as a good and convenient marker of quantum 
phase transition. In a strongly correlated system, the 
reduced density matrix of a local block involves the bi- 
partite correlation in the system as the consequence of 
the superposition principle in quantum mechanics. At 
the quantum critical point, quantum fluctuations extend 
to all scales, thus it may contribute a singular part in 
the local entanglement. The idea of using the local en- 
tanglement to locate transition point is especially useful 
in numerical calculations [7J, usually the data obtained 
from a relative small size lattice is sufficient to detect the 
quantum phase transitions. 

In this paper, we investigate the scaling behavior of 
the two-site local entanglement near the quantum phase 
transition point in ID transverse field Ising model. Since 
the transition type of the Ising model is different from 
that in ID XXZ model, the entanglement shows quite 
different behavior, such as scaling and singularity. We 
calculate the entanglement and its first order derivative 
with respect to the coupling constant numerically. Based 
on the exact solution [la . Il7| , we also obtain explicit ex- 
pressions which clearly exhibit the logarithmic divergence 
of local entanglement's derivative at the critical point. 
Therefore, by studying the two-site local entanglement, 
we hope to have a deep understanding on the another 
type of quantum phase transition, as represented by ID 
Ising model. 
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The Hamiltonian of ID transverse field Ising model 



reads 
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where ai(a x ,a v , a z ) are Pauli matrices at site i, A is an 
Ising coupling in unit of the transverse field, and the 
periodic boundary conditions are assumed. In the ba- 
sis spanned by the eigenstates of {<r z }, the Hamiltonian 
changes the number of down spins by two due to the 
term afcrf +1 in its expression, so the whole space of sys- 
tem can be divided by the parity of the number of down 
spins. That is the Hamiltonian and the parity operator 
P = Y[j Oj can be simultaneously diagonalized and the 
eigenvalues of P is ±1. Then it can be proved that for 
a finite system, the ground state in the whole region of 
A > is non-degenerate. 

We confine our interest to the entanglement between 
two neighboring sites and rest of the system, so we need 
to consider the reduced density matrix of two local sites. 
According to the parity conservation, the reduced density 
matrix of two spins in this system takes the form 
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in the basis {| TT>, | tl), | IT), | II)}- The elements in the 
density matrix p^j can be calculated from the correlation 
functions 
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In a finite system, the ground state is a pure state, so the 
von Neumann entropy E v , i.e. 



E v = -tr [pijlog^pij)] 



(4) 



measures the entanglement between the sites i, j and the 
rest N — 2 sites of the system. 

The correlation functions in Eqs ([3]) can be calculated 
exactly. The mean magnetization in the ground state is 
given by 16] 



{<>*) 



1 ^ (1 — Acos0) 



(5) 



where lo^ is the dispersion relation, 

U(l> = J~l + A 2 - 2Acos((/)g), 
(f> q = 2nq/N, 



(6) 
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FIG. 1: (color online) The two-site local entanglement as a 
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FIG. 2: (color online) The first order derivative of the two-site 
local entanglement as a function of A for various system size. 



where q is integer (half-odd integer) for parity P = 
— 1(+1). The two-point correlation functions are calcu- 
lated as [17 1 



T y~y\ — 



1 ^-^ A — cos 
N 1^ 
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A cos(2</>) — cos ( 



(7) 



Therefore, we can calculate the two-site local entangle- 
ment for arbitrary finite system directly. 

As is well known, in a ring with infinite sites, the phase 
diagram of the transverse field Ising model is divided into 
two phases by a transition point A = 1. For A < 1, the 
ground state of the system is a paramagnet. Its phys- 
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FIG. 4: The scaling behavior of the maximum point of 
dE v /d\. The circle line is obtained by numerical derivative, 
and the square line by analytical expression directly. 

ical picture can be understood from the case of A = 
at which all spins are fully polarized along z direction. 



While for A > 1, the strong Ising coupling introduces 
magnetic long-range order of the order parameter (cr x ) 
to the ground state. Therefore, the transition in this 
model is of type from order-to-disorder, and is quite dif- 
ferent from the transition in some other models. In the 
ID XXZ model, it has been shown that the two-site local 
entanglement reaches a local maximum at the transition 
point A = 1 which witnesses a transition of type from 
order-to-order, and becomes singular at another transi- 
tion point A = — 1 due to the ground-state level-crossing. 
However, in the ID Ising model, no similar properties can 
be observed directly from the entanglement, as shown in 
Fig. [TJ As suggested by the classification scheme for 
the quantum phase transition based on the pairwise en- 
tanglement [1« ] , we take the first order derivative of the 
entanglement with respect to the coupling A. The results 
are shown in Fig. [2] We are happy to notice that the 
first order derivative around the critical point becomes 
sharper as the system size increases, and is expected to 
be divergent in an infinite system. 

In the critical phenomena, the most important themes 
are the scaling and universality. In a finite sample, a 
physical quantity is a smoothly continuous function if 
there is no ground-state level-crossing. However, it is ex- 
pected that the anomalies will become clearer and clearer 
as the size of sample increases. The relevant study is the 
finite size scaling. For this system, in order to study 
the scaling behavior the maximum point of dE v /dX, we 
locate the maximum point for a given system size numer- 
ically and define the corresponding location as A m . From 
Fig. [2l though there is no divergence when N is finite, 
the anomalies are obvious. The position of the maximum 
point A m scales as A c — A TO oc 1 /TV 3 / 2 . 

The results of dE v /d\ at A m up to 2000 sites are shown 
in Fig. 31 Though we can almost conclude logarithmic 
divergence of dE v /dX at the transition point from the 
figure, it is still very difficult to take numerical derivative 
for a very large system, such as N = 10 s . In order to 
make it more confirmative, we now give a rigorous prove 
that dE v /dX obeys a logarithmic behavior at the critical 
point. 

In order to calculate the entanglement, we first need to 
calculate the eigenvalues of the reduced density matrix 
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Then the first order derivative of the entanglement takes At the critical point, the eigenvalues Ci converges quickly 

as the system size increases 



dE, 
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Then after some complicated calculations, one can obtain 

(11) 
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The main contribution to the above expression in the 
large N limit arise from the summation around zero point 
of 6. Because 
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then we obtain, in the large N limit 
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where Af is a constant 



A? 



J . 3tt 2 + 4tt - 16 



" 2tt l0g2 3tt 2 - 4tt - 16 



3 , 3tt 2 

^=— logo 

2VT3tt 2 3tt 2 



16 



4%/13tt + 16 



(14) 



as has been shown in Fig. 2J On the other hand, for 
an infinite system, dE v /d\ diverges on approaching the 
critical point as 



dE v 
~~dX 



= A% In | A - A c 



const, 



(15) 



as has been shown in Fig. [5] Here A® is not completely 
independent of A, and Af = Af at A = 1. 

According to the scaling ansatz the two-site en- 
tanglement, considered as a function of the system size 
and the coupling, is a function of N 1 / V {X — A m ). In the 
case of logarithmic divergence, it behaves as dE v /dX — 
dE v /dX\ x =x m ~ Q[N^ V (X- A m )], where Q{x) oc hxx 
for large x. In Fig. [SI we perform the scaling anal- 
ysis and find the entanglement can be approximately 
collapse to a single curve for v ~ 0.98. This be- 
havior is quite different from that of concurrence. (20j 
which diverges as dC/d\\\=\ m oc — 0.27021nA^ and 
dC/dX\\ = \ c oc (8/37T 2 ) In |A — A c | respectively. Because 
of 8/37T 2 ~ 0.2702, the ratio of constants in the In term 
is unit. This fact leads to v = 1 for the concurrence. [20j 
Actually, though Osterloh et al only gave a numerical 
constant 0.2702 for the logarithmical divergence of the 
concurrence with increasing system size, it also can be 
exactly obtained that the constant is just 8/37T 2 with the 
above procedure. However, for the local entanglement, 
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FIG. 5: Logarithmical divergence of dE v /d\ around the crit- 
ical point of an infinite system. 
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FIG. 6: (color online) The finite-size scaling analysis for the 
case of logarithmic divergence. The local entanglement, con- 
sidered as a function of system size and the coupling, collapse 
on a single curve for various system size. 



since Af slightly depends on A, we have v ~ 0.98 ac- 
cording to numerical analysis. 

In summary, we have investigated the scaling behavior 
of the two-site local entanglement in a quantum phase 
transition of ID transverse field Ising model. Different 
from the entanglement in the transition of ID XXZ model 
at A = 1 which is of type from order-to-order, the en- 
tanglement in ID transverse field Ising model is not a 
maximum at the critical point where the quantum phase 
transition is of the type from order-to-disorder. However, 
its first order derivative with respect to the coupling be- 
comes singular around the critical point as the system 
size increases. We show the logarithmic divergence both 
numerically and analytically. 

Note added: During the preparation of the manuscript, 
we notice that a work on the two-site entanglement in 
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